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Abstract 

The mechanisms of cold- and pressure-denaturation of proteins are matter of debate and are 
commonly understood as due to water-mediated interactions. Here we study several cases of 
proteins, with or without a unique native state, with or without hydrophilic residues, by means 
of a coarse-grain protein model in explicit solvent. We show, using Monte Carlo simulations, that 
taking into account how water at the protein interface changes its hydrogen bond properties and 
its density fluctuations is enough to predict protein stability regions with elliptic shapes in the 
temperature-pressure plane, consistent with previous theories. Our results clearly identify the 
different mechanisms with which water participates to denaturation and open the perspective to 
develop advanced computational design tools for protein engineering. 
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Water plays an essential role in driving the folding of a protein and in stabilizing the 
tertiary protein structure in its native state [l|, 2|. Proteins can denaturate—unfolding 
their structure and loosing their activity—as a consequence of changes in the environmental 
conditions. Experimental data show that for many proteins the native folded state is stable 
in a limited range of temperatures T and pressures P JsHsI and that partial folding is T- 
modulated also in “intrinsically disordered proteins” |9|. By hypothesizing that proteins 
have only two different states, folded (/) and unfolded (nh and that the / <—> u process 
is reversible at any moment, Hawley proposed a theory [1^ that predicts a close stability 
region (SR) with an elliptic shape in the T — P plane, consistent with the experimental data 


111 . 


Cold- and P-denaturation of proteins have been related to the equilibrium properties of 
the hydrat 


versial 


ion water 


12 


24 


23|. However, the interpretations of the mechanism is still contro- 


331]. High-T denaturation is easily understood in terms of thermal fluctuations 


that disrupt the compact protein conformation: the open protein structure increases the en¬ 
tropy S minimizing the global Gibbs free energy G = H — TS, where H is the total enthalpy. 
High-P unfolding can be explained by the loss of internal cavities in the folded states of pro¬ 
teins , while denaturation at negative P has been experimentally observed 3^ and simu¬ 
lated [3^ recently. Cold- and P-unfolding can be thermodynamically justified assuming 
an enthalpic gain of the solvent upon denaturation process, without specifying the origin of 
this gain from molecular interactions 3q . Here, we propose a molecular-interactions model 


for proteins solvated by explicit water, based on the “many-body” water model 


29 


37 




We demonstrate how the cold- and P-denaturation mechanisms can emerge as a competition 
between different free energy contributions coming from water, one from hydration water 
and another from bulk water. Moreover, we show how changes in the protein sequence af¬ 
fect the hydration water properties and, in turn, the stability of the protein folded state—a 
relevant information in protein design 25]. 

The many-body water model adopts a coarse-grain (CG) representation of the water 
coordinates by partitioning the available volume V into a fixed number Nq of cells, each with 
volume V = V/No > vq, where vq is the water excluded volume. Each cell accommodates 
at most one molecule with the average 0-0 distance between next neighbor (n.n.) water 
molecules given by r = To each cell we associate a variable rij = 1 if the cell i is occupied 
by a water molecule and has Vq/v > 0.5, and n* = 0 otherwise. Hence, rii is a discretized 
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density field replacing the water translational degrees of freedom. The Hamiltonian of the 
bulk water 

U{r,j) - - J.iVeoop (1) 

ij 


has a hrst term, summed over all the water molecules i and j at 0-0 distance Vij, accounting 
for the van der Waals interaction, with U{r) = oo for r < tq = = 2.9 A (water van der 

Waals diameter), U{r) = 4e[(ro/r)^^ — (ro/r)®] for r > tq with e = 5.8 kJ/mol and U{r) = 0 
for r > rc = 6ro (cutoff). 

The second term represents the directional component of the hydrogen bond (HB), with 
J/4e = 0.3 42|, iV^g = X](p) number of bulk HBs, with the sum over n.n., where 

(Tij = 1,... ,g is the bonding index of molecule i to the n.n. molecule j, with 5ab = 1 if 
a = 6, 0 otherwise. Each water molecule can form up to four HBs that break if riirij = 0, i.e. 
Tij > 21/3^0 = 3.6A, or OOH > 30°. Hence only 1/6 of the entire range of values [0, 360°] for 
the OOH angle is associated to a bonded state. Thus we choose g = 6 to account correctly 
for the entropy variation due to HB formation and breaking. 

The third term, with N^oop = k)i^(^ik,o-u^ where {l,k)i indicates each of the six 

different pairs of the four indices cr^ of a molecule i, accounts for the HB cooperativity due 
to the quantum many-body interaction 43| and leads to the low-P tetrahedral structure 
44l |. We choose Jo-/4e = 0.05 <C J, to guarantee an asymmetry between the two HB terms. 

Increasing P partially disrupts the open structure of the HB network and reduces v toward 
uq. We account for this with an average enthalpy increase Pn^g per HB, where 4 b/^o = 0.5 
is the average volume increase between high-p ices VI and VHI and low-p (tetrahedral) ice 
Ih. Hence, the total bulk volume is = Nvq + fVggnjj’AWe assume that the HBs do not 
affect the n.n. distance r, consistent with experiments |4^, hence do not affect the U{r) 
term. 

Next we account for the effects of protein-water interaction. Experiments and simulations 
show that near a hydrophobic (<h) residue water-water HBs are more stable then in bulk 
[l5, 45-47] with stronger water-water correlation j^|. We model this by replacing J of 
Eq. ([I]) with > J for HBs at the 4) interface. This choice, accor ding to Muller c 

ensures the water enthalpy compensation during the cold-denaturation 49]. 


The interaction at the 4> interface affects the hydration wa 


Some works suggest a decrease of interfacial water density 


50 


nr density and fluctuations. 


53| . while recent simulations 
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show an increase of density in the hrst hydration shell of any solute 


compressibility near $ solutes with size > 0.5 nm for water j26l . l48l . 


54j and an increase of 


55l | or water-like solvents 


56l | with respect to bulk. Increasing P induces a further increase of 


density and reduces the 


effect of the $ interface on the compressibility of the hydration shell 


26 


48j . We incorporate 


this behavior into the model by using the following considerations. From the equilibrium 
condition for the thermodynamic potential of hydration water and the coexisting vapor at the 
$ interface at hxed T, according to the Eq. (2) of Ref. 57|, we deduce ~ (P—P*)“\ 

where is the volume per hydration water molecule and P* < 0 is the equilibrium vapor 
pressure at the given T. If we attribute this P-dependence to the interfacial HB properties 


,(‘f) 

water-water HB at the <F interface is 


■^HB ~ ~ "^o) and expand it as a power series in P, the average volume change per 


4bMb,o = l-hP + k2P^ - hP^ + 0(P") (2) 

where volume change associated to the HB formation in the $ hydration shell 

at P = 0, fcj > 0 Vi and limp^oo'y^ffi = 0. Hence, the total volume V is 

y = y(b) + = y(b) + (3) 

where and are the <F hydration shell volume and number of HBs, respectively. 
Because we are interested to small values of P, i.e. near the biologically relevant atmo¬ 


spheric pressure, we include in our calculations only the linear term in Eq.(|2]) [58|]. In the 
following we hx ki = lno/4e, = 0.5 and J<^/J = 1.83. Our results have 

minor qualitative differences by including up to the third order in Eq. ([2]) or by changing up 
to 50% the parameters. 


Because our goal here is to calculate the water contribution to dena 
the protein as a self-avoiding $ homopolymer, without internal cavities 


uration, we model 


59l |. whose residues 


occupy n.n. cells with no resic 
CG approaches to the problem 


ue-residue interaction but the excluded volume, as in other 
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m 


60l | . This implies that the protein has several “native’' 


states, all with the same maximum number Umax of residue-residue contacts. To simplify the 
discussion, we initially neglect energetic contributions of the water-<F residue interaction. 



visualization and understanding of our results. Comparisons with our preliminary results in 
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Temperature T 


FIG. 1. P — T stability region (SR) of the protein from MC simulations. Symbols mark state 
points with the same average residue-residue contact’s number rerr/umax = 30%, 40%, 50% and 
70%. Elliptic lines are guides for the eyes. The “glass transition” line defines the temperatures 
below which the system does not equilibrate. The spinodal line marks the stability limit of the 
liquid phase at high P with respect to the gas at low P; ks is the Boltzmann constant. 


3D do not show qualitative changes, mainly because the number of n.n. water molecules is 
four both in 2D and 3D for the tendency of water to form tetrahedral structures in 3D. 


We consider that the protein is folded if the average number of residue-residue contacts 
Wr > 50% Umy. We hnd an elliptic SR (Fig{T]), consistent with experiments and the Hawley 
theory 10, lll|, with heat-, cold-, and P-unfolding. The elliptic shape is preserved when 
we change the threshold for showing that the / <—)■ n is a continuous process. In the 
SR the folded protein (Fig. [2^) minimizes the number of hydrated <F residues, reducing the 
energy cost of the interface, as expected. 


First, we observe that the model reproduces the expected entropy-driven f <—> u for 
increasing T at constant P (Fig. [2b). The entropy S increases both for the opening of the 
protein and for the larger decrease of HBs. 


Upon isobaric decrease of T the internal energy dominates the system Gibbs free energy 
(Fig. [2b). However, saturates at T lower than the SR, therefore the only way for 
the system to further minimize the internal energy is to increase i.e. to unfold the 

protein. Hence the cold denaturation is an energy-driven process toward a protein state that 
is stabilized by the increased number of HBs in the hydration shell. 
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FIG. 2. Typical configurations of a hydrated protein made of 30 residues (in green): (a) folded 
at the state point (Ffc^/de, Puo/4e) = (0.25,0.1) and unfolded (b) at high-T (0.9,0.1); (c) at 
low-T (0.1,0.1); (d) at high-P (0.25,0.6); (e) at low-P (0.25,—0.3). Left panels: Water molecules 
with/without HBs are represented in blue/white and bulk/interfacial HBs in blue/red. Right 
panels: Color coded water density field (from black for lower p to yellow for higher p) calculated 
as vqp^^'^ = uo/(uo + where A = b, <I>,gand n^g ^ is the number of HBs associated to the 

water molecule z, with rigg ^ = IVgg. 











































Upon isothermal increase of P the protein denaturates (Fig. [2]i). We hnd that this change 
is associated to a decrease of and a small increase of that lead to a net decrease 
of V at high P, as a consequence of Eqs. ([2]) and ([2D, and an increase of internal energy. At 
high P the PU-decrease associated to the / —)■ u process at constant T dominates over the 
concomitant internal energy increase, determining a lower Gibbs free energy for the u state. 
Hence water contribution to the high-P denaturation is density-driven, as emphasized by 
the increase of local density near the unfolded protein. 

Finally, upon isothermal decrease of P toward negative values (Fig. [2^), the enthalpy 
decreases when the contribution (Pn^^ — </$)A'hb decreases, i.e. when increases. 

Therefore we hnd that under depressurization the denaturation process is enthalpy-driven. 

From the Clapeyron relation dP/dT = AS/AV applied to the SR we expect that 
the / i—> u process is isochoric at the SR turning points where dT/dP\sB. = 0, while is 
isoentropic at the turning points where (9 P/c1T|sr = 0. In particular, at any T and P the 
volume change in the / —> u process is given by 


AV' = 14 - ~ - A^iPjAAf®. 


(4) 


We estimate the Eq. (jlD calculating the average volume W and Vf in a wide range of T 
and P, equilibrating water around a completely unfolded protein state and a completely 
folded state (with n„ = Umax)- Consistently with the Hawley’s theory 1^, we hnd that 
the T-denaturation is accompanied by a positive entropy variation AS* > 0 at high T 
and an entropic penalty AS < 0 at low T, while the P-denaturation by a decrease of 
volume AV < 0 at high P and an increase of volume AV > 0 at low P (Fig. [3D. In 
particular, at P = 0.3(4e/no), corresponding to ~ 500 MPa, we hnd AV ~ —2.5no, hence 
|PA1/| = 0.75(4e) ~ 17 kJ/mol, very close to the typical reported value of 15 kJ/mol Q • By 
varying the parameters and J<j, we hnd that the hrst is relevant for the P-denaturation, 
as expected because it dominates Eq. ([HD, while the second ahects the stability range in 
T. Both combine in a non-trivial way to regulate the low-T entropic penalty. We test our 
results including a small water-<h residue attraction and hnd no qualitative diherences but 
a small change in the T-range of stability of the folded protein. 


N 


ext, we study the case of a protein with hydrophobic (<F) and hydrophilic (C) residues 
61| . with a residue-residue interaction matrix Aij = e^r if residues i and j are n.n. 
in the unique native state, 0 otherwise. Water molecules interact with energy < J 
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FIG. 3. Volume change AV for the / —)■ u process in the T — P plane. Color coded AV (black 
for negative, yellow for positive) is in vq units. Solid lines connect state points with constant AV. 
Black points mark the SR. The locus AV = 0 has a positive slope and intersects the SR at the 
turning points with dT/dP\sR = 0. The dashed line, connecting the points with dP/dT\sR = 0, 
corresponds to the locus where AS = 0 and separates state points with AS" > 0 (high T) from 
those with AS" < 0 (low T) at the / —> u process. The white symbol marks the error in the 
dashed-line slope estimate. 



Temperature T [48/^^] 


FIG. 4. The SR for the heteropolymer with a unique native state. We set e„/J = 0.7, e*,# = 0, 
— 1-17, J$/J = 1.3, J(^/J = 0.5, Uyg = 0, with all the other parameters as in Fig. [TJ We 
test that changing the parameters, within physical ranges, modifies the SR, reproducing a variety 


of experimental SRs 


11[, but preserves the elliptic shape. 

















and ew,c > J with n.n. $ and C residnes respectively, accounting for the polarization of 
the solvent near the C, residues. The polar ^ residues disrupt the tetrahedral order of the 
surrounding water molecules. Thus we assume that a ^ residue j and a n.n. water molecules 
i form a HB when the latter has (jjj in the state q^p = 1,..., g preassigned to j. Finally, 
we consider that water-water enthalpy in the hydration shell is Hx x = —J\ + both 

molecules are n.n. to the same type of residue or Hx^^j. = (Hx x + if the n.n. residues 

belong to different types, with A,/i = and J^< J (62| (FigJl]). 

Despite the complexity of the heteropolymer model, we find results that are similar 
to the homopolymer case, with the qualitative difference that with (^-residues we hnd a 
locus AD = 0 with negative slope and increased stability toward (i) cold- and (ii) P- 
denaturation. In particular, for our specific choice of parameters, for the heteropolymer the 
cold denaturation at P = 0 occurs below the glass transition, instead of ~ 50% above as for 
the homopolymer. Furthermore, the SR against P is ~ 2 times larger with (^-residues than 
without. This comparison suggests that the water contribution is relevant to the / < —)■ u 
independently on the residue sequence, although the residue-residue interactions increase 
the stability of the folded state. 

In conclusion, our model for protein folding reproduces the entire protein SR in explicit 
solvent and allows us to identify how water contributes to the T- and P-denaturation pro¬ 
cesses. The model is thermodynamically consistent with Hawleys theory but, in addition, 
allows for intermediate states for the / <—)■ u process. We find that cold denaturation is 
energy-driven, while unfolding by pressurization and depressurization, in addition to other 
suggested mechanisms 33(], are density- and enthalpy-driven by water, respectively. For 
these mechanisms is essential to take into account how the protein-water interactions af¬ 
fect the stability of the water-water HB and the water density in the hydration shell. In 
particular, both properties control the low-T entropic penalty. Our results are qualitatively 
robust against modihcation of the model parameters, within physical ranges, and the model 
is computationally efficient thanks to the adoption of a CG water model, representing a step 
towards the development of a theoretical and computational approach for protein design 
and engineering. 
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